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The cytoskeleton is an important subsystem of cells that is involved for example in cell division 
and locomotion. It consists of filaments that are cross-linked by molecular motors that can induce 
relative sliding between filaments and generate stresses in the network. In order to study the effects 
of fluctuations on the dynamics of such a system we introduce here a new class of driven diffusive 
systems mimicking the dynamics of active filament bundles where the filaments are aligned with 
respect to a common axis. After introducing the model class we first analyze an exactly solvable 
case and find condensation. For the general case we perform a mean-field analysis and study the 
behavior on large length scales by coarse-graining. We determine conditions for condensation and 
establish a relation between the hopping rates and the tension generated in the bundle. 
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I. INTRODUCTION 

> , 

' Our current understanding of non-equilibrium systems owes a lot to the studies of driven diffusive systems 0. 
These are models of particles hopping from site to site on a lattice where the hopping rates do not fulfill the condition 
of detailed balance. Even though the rules governing the stochastic dynamics are often very simple, these models 

S show a variety of phenomena that are unknown to systems at thermodynamic equilibrium. In particular, even in one 
spatial dimension they can exhibit phase separation and phase transitions. 
"-^J , Biology provides ample examples of systems that are maintained out of thermodynamic equilibrium and physical 

■ studies of these systems are essential for understanding their behavior. In contrast to the driven diffusive systems 
mentioned before, they in general do not follow simple rules. This can be illustrated by the cytoskeleton which is a 

i network of filamentous proteins, mainly actin-filaments and microtubules 0,0,0]. The cytoskeleton plays an essential 

role in many vital cellular processes as, for example, cell division and locomotion. Cytoskeletal filaments are polar 

■ objects, because they have two structurally distinguishable ends. Other proteins associated with the cytoskeleton 
control filament assembly and disassembly. Furthermore, they can form cross-links between filaments. In addition 
to passive cross-links there are active cross-links formed by motor proteins. These are molecules that can move 
directionally along filaments. Their motion is driven by the hydrolysis of Adenosine- Tri-Phosphate during which a 

^ phosphate group is cleaved from the molecule. The chemical energy thus liberated is then transformed by the motor 
proteins into mechanical work. The direction of their motion is determined by the orientation of the filament. If a 
motor protein connects two filaments its motion results in relative motion between the filaments and creates stresses 
. in the filament network. 

Now each of the proteins involved in the dynamics of the cytoskeleton is already a complicated system on its own, 
consisting of many atoms. However, the separation of time scales allows phenomenological descriptions of the slow 
modes. Essential features of cytoskeletal dynamics can be studied using simple models. For example, driven diffusive 
systems have been used to study the growth of filaments against external forces the motion of single motors on a 
filament 0, Q , as well as collective effects of motor molecules interacting through excluded volume effects 0, 0, . 

Motivated by the dynamics of bundles of polar filaments interacting with molecular motors, we study in this work 
a new class of driven diffusive systems. In contrast to commonly studied models, elementary events consist here of a 
correlated motion of two particles reflecting the effects of filament motor interactions. After defining the model class 
we study a sub-class that allows for an exact solution of the steady state. We identify criteria for condensation in 
these models. In a filament bundle, condensation corresponds to bundle shortening. We then analyze cases without 
JH \ exact solution by employing a mean-field approximation. Finally, we look at a coarse-grained version and relate the 
model to phenomenological descriptions of active filament bundles. In this way we can connect the kinetic hopping 
rates defining the driven diffusive system to the stress generated in a filament bundle. 



o 
in 



-a 



'Electronic address: pk.mohanty@saha.ac. in 
^Electronic address: karsten@pks.mpg.de 



2 



II. ACTIVE FILAMENT BUNDLES 

We will now briefly review a microscopic physical model of the filament dynamics in a bundle induced by motor 
molecules. Motivated by this model we will then define the processes that determine the class of driven diffusive 
systems to be studied below. 

A. The minimal model 

Consider a bundle of polar filaments in which all filaments are aligned with respect to a common axis. Each filament 
has one of two possible orientations which in the following will be referred to as "plus" and "minus" , respectively. 
Filament pairs are set into relative motion by the action of molecular motors which form active cross-links. As a 
consequence, the bundle can shorten and filaments can be separated according to their polarity [TTL Il2| . 

For sufficiently large filament densities, the state of the bundle can be represented by two densities corresponding 
to the distributions of filaments oriented in one way and the other. They are denoted c + and c~ , respectively. Due 
to mass conservation, the time evolution of the densities is governed by continuity equations 

dtc* + aj* = o (i) 

The currents reflect various active processes in the bundle, for example the action of motors connecting filament 
pairs. If filaments are created and destroyed or change their orientation, then source terms have to be added on the 
right hand side of Eq. QJ. In the following, these processes will be neglected such that the numbers of plus- and 
minus-filaments are both conserved independently. 

Different approaches have been developed to derive expressions for the filament currents. Phenomcnological equa- 
tions have been derived based on systematic expansions of the currents in terms of the filament densities and their 
derivatives respecting the symmetries of the system 0, \l$L lri| . The same approach has been used for de- 

scribing the hydrodynamics of actively propelled particles like bacteria which share a number of essential features 
with active filament systems |2jj. The virtue of these approaches is that they are independent of many mi- 
croscopic details of the active processes, most of which are only partially known at present. In more microscopic 
descriptions, expressions for the currents follow from an analysis of the action of active cross-links formed by motor 
molecules |1 |H |H H| . 

We will now discuss further the approach presented in Ref. |22| as it will form the basis for the driven diffusive 
systems discussed later. It is based on an analysis of the momentum flux in the bundle. Since the motion of filaments 
in a bundle typically occurs with velocities in the order of /im/min, the dynamics is dominated by friction and inertial 
terms can be neglected. The conservation of momentum is thus expressed by a balance of forces. Summing up all 
the forces acting on a filament connected to other filaments by active motors and moving against an immobile fluid 
leads to a general expression for the currents in terms of the internal forces. Within the framework developed in 
Ref. |22j it is assumed that the dominant contributions to the currents result from filament pair-interactions, which 
is for example valid in the case of a low motor density. In the simplest case, the motor density is not considered as a 
dynamic variable, such that it only enters in the parameters describing the dynamics. 

For a simple model of the internal forces between filament pairs, the current of plus-filaments is given by |2Hl22| 

j+(x) = -ad x / d^[c+(x + 0-c+(x-0}c + (x) 



+P8 X / d£ c~(x + - Dd x c+{x) + v tl .c+{x) (2) 

J -I 

This current defines the minimal model of active filament bundles. The first term describes the contribution to 
the current resulting from interactions between filaments of the same orientation, while the second term accounts for 
interactions between filaments of opposite orientation. The corresponding coupling constants are a and f3, respectively. 
The integrals reflect that filaments can interact whenever they overlap, where £ is the filament length which is taken 
to be the same for all filaments. The terms resulting from the action of active cross-links have a simple interpretation. 
The interaction of filaments of the same orientation is such that, for a > 0, a filament pair tends to increase its 
overlap, while for filaments of opposite orientation the distance between their plus ends decreases if (3 > 0, see Fig. ^ 

The following term is a diffusion current that captures effects of fluctuations in the system. Since the fluctuations 
are not only thermal, the effective diffusion constant D does not satisfy an Einstein-relation. Finally, the convective 
part v tI c + of the current has been introduced in Ref. |25| to describe filament treadmilling. In this process filament 
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FIG. 1: The three processes governing the particle dynamics, a) Particles of the same kind separated by one site, hop 
simultaneously to the middle site with rate u. On the right the corresponding process of two filaments with the same orientation 
aligning their plus-ends is illustrated; the coupling parameter in the minimal model is a. b) The opposite event of two particles 
of the same kind leaving a site to different neighboring sites occurs with rate v. c) Two neighboring particles of different kinds 
exchange their position with rates w^ 1 . On the right the corresponding process of sliding two filaments of opposite orientation 
is illustrated; the corresponding coupling parameter in the minimal model is (5. d) Particles hop to the left or right neighboring 
sites with rates cLl r- 



subunits are added at the plus-end and removed with the same rate at the minus-end. This leads to an apparent 
motion of the filament in the direction of the plus-end while the length of the filament remains constant. 

The expression for the current of minus-filaments is obtained from expression J2J by interchanging the superscripts 
+ and — and by replacing (3 with — (3 as well as Wtr by — v tv . This assures that the currents transform correctly under 
space inversion which also transforms plus filaments into minus filaments and vice versa. Furthermore, the currents 
respect the conservation of momentum. 

An analysis of the dynamic equations reveals that the important parameter for the dynamics is a, the coupling 
between filaments of the same orientation. If it exceeds a critical value, a homogeneous distribution of filaments 
becomes unstable [^[i^. Provided that (3 = 0, filaments then aggregate at some position. If on the contrary /3 ^ 0, 
then oscillatory solutions in the form of traveling waves can be generated [26]. The critical value a c is inversely 
proportional to the filament density Cq of the homogeneous state, a c oc Cq . 

We will now use these elementary processes as a motivation for defining a class of driven diffusive systems capturing 
essential parts of the dynamics of active filament bundles. 



B. Stochastic dynamics of interacting particles 

Consider a one-dimensional lattice of L sites. On the lattice there are two kinds of particles, referred to as "plus" 
and "minus", which represent the two possible orientations of the filaments in a bundle. Each site i with i = 1, . . . , L 
accommodates nf + > particles, where nf are the occupation numbers of plus- and minus-particles. The 
dynamics of the particles is determined by processes corresponding to the different terms appearing in the current 
©, see Fig.Hl 

First of all, two particles of the same kind located respectively on sites i — 1 and i + 1 can simultaneously hop to site 
i. This process corresponds to the dynamics of a pair of equally oriented filaments in a bundle that are cross- linked 
by a motor. The action of the motor tends to increase the overlap of the filaments, see Fig. right. The rate 
of such an event depends on the number of particles on sites i — 1 and i + 1 and is denoted by u(nf_ 1 , nf +1 ) with 
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u(n, m) = u(m, n). Note, that the function u is the same for both kinds of particles 1 . For the following analysis, it will 
turn out to be convenient to allow also for the opposite events to occur even though there is no direct analog in the 
minimal model: Two particles on site i can separate and move to sites i — 1 and i + 1, respectively. The corresponding 
rate v depends on the number of particles at i, i.e. u(rij). 

The next process reflects the action of a motor on a pair of oppositely oriented filaments, illustrated in Fig. ^i, 
right. A plus particle at i will exchange its position with a minus particle at i + 1 with a rate w + that depends on 
nf and n~ +1 , i.e. w + (nf,n~ +1 ). The opposite process occurs with rate w~(n~,nf +1 ). Motivated by the behavior 
of the term proportional to /3 in Eq. @ under the combined operation + «-> — and space inversion, we require 
(m, n) — w T (n, m) . 

Furthermore, we allow for processes similar to diffusion: particles can spontaneously hop to the left or right 
neighboring site with rates df^ and <i^, respectively, each of which can depend on the number of particles at the site. 
Motivated by the behavior of the filament system under space inversion, we require d\ — and d^ = d^. Note, that 
if d~^ ^ d~^ then there is a spontaneous flux of particles similar to the convective term uti-c + in Eq. Pfl. Finally, the 
rates have to respect that a process can only take place if there is a particle at the donating site(s), i.e., u(n, n') = 
if either n = or ri = 0, u(0) = u(l) = 0, and df R (0) = 0. 

If only one kind of particles is p resent and if u = v = 0, the model reduces to a class of hopping models known 
as Zero Range Process (ZRP) [22|. In these models, the rate of a particle leaving a site depends on its occupation 
number but is otherwise unconstrained. In spite of this simplicity, phase transitions are found in this class even in one 
dimension. The ZRP allows for an analytical solution of the steady state. This property is intimately linked to the 
factorization of the steady state probability distribution P({rii}) of having a certain configuration with rn particles 
at site i = 1, . . . , L. Explicitly, 

L 

P({n i }) = Z- l Hf(n j ), (3) 

3=1 

where f(n) depends on the hopping rate. The factor Z assures normalization of the probability distribution, depends 
on the length of the lattice L and the number of particles N in the system. For example, if di,{n) = pd(n) and 
dji(n) — qd(n) with arbitrary p > and q > 0, then 

W - A ^Ty (4) 

m=l y 1 

If the rate d(n) goes to zero asm oo, the system shows condensation, i.e., in the steady state there will be a finite 
number of sites with a macroscopic occupation number [2^ ■ 

We will now analyze the case of non-vanishing u and v, where we focus on the case of one kind of particles. First, 
we will study a case that allows for an analytic solution of the steady state probability distribution and shows that 
the the model can show condensation. Then we discuss the general case by using a mean-field approximation. 



III. FACTORIZED STEADY STATE 



Consider the case of N particles of one kind that move on a periodic lattice of L sites. The hopping rates are 
specified by 

u(m,n) = rw(m)w(n) (5) 

v{n) = rw(n)w(n — 1) (6) 

dh{n) = qw(n) (7) 

dfl(n) = pw(n) (8) 

where w is an arbitrary function of n with w(0) = and p, q, and r are positive parameters. As can easily be checked, 
these rates fall into the class introduced in the previous section. For this process, the steady state probability 



1 This process is similar to the one introduced in Ref. |2 ll where the "interaction range" was larger and u was chosen to be constant. 
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distribution factorizes and has the form 10. Indeed, the steady state of the probability distribution P is determined 
by 

L 

^2 [2u(nj_i + 1, n i+1 + 1)P(. . . , n^i + 1, m - 2, n i+1 + 1, . . .) 

-{u(rii_ 2 ,«i) + u(ni,n l+2 ) + 2v(n l )}P(. . . , n;_i,ni, n i+ i, . 

+w(ri i _i + 2)P(. . . , rii_ 2 - 1, + 2, - 1, . 
+v(n i+1 + 2)P(. . . , 7ij - 1, + 2, n i+2 - 1 . 
+d R (n i - 1 + 1)P(. . .,n,i-i + l,nj - 1, . . .) - d L (ni)P{. . .,m-i,rH, . 
+d L (n i+1 + 1)P(. ..,m-l, n i+1 + 1, . . .) - d R {ni)P{. . . ,n u n i+1 , . . .)] = 0. 

(9) 

This condition is clearly fulfilled if each process occurs with the same rate as the corresponding opposite process. 
Inserting the ansatz this detailed balance yields 3L conditions 



u(m-i + 1, n i+ i + l)f(m-x + l)/(rij - 2)/(n i+ i + 1) = 
dh{ni + i + l)f(m+i + l)f(rii - 1) = 
d R {rii-i + l)/(nj_i + l)f(rii - 1) = 

with i = 1, . . . , L. For the special choices ©-(JSl), they are solved by 

. f(n + 1) 
win + 1) — t^— r — — const. 

f( n ) 

Without loss of generality we choose the constant to be 1 and find 

n 1 
/(n)=II 



v(ni)f{ni^i)f(ni)f(n i+ i) 
dh{ni) j '(rii_i) j '(m) 
dR(ni)f(m)f(n t+1 ) 



1 w(n) 



(10) 



(11) 



(12) 



which formally is the same expression as for the ZRP. 

The similarity of the steady state probability distribution with that of the ZRP allows us to immediately carry over a 
number of results to our system. Of particular interest in the present context are the conditions for which the particles 
condensate. Condensation means that in the "thermodynamic limit" , where N, L — > oo with p — N/L — constant, a 
single site is occupied by macroscopic number of particles. Three cases can be distinguished as follows : 

1. w(n) — > oo as n — * oo: No condensation, the distribution of particles will be homogeneous in the limit N — > oo. 

2. w(n) — > as n — > oo: Condensation occurs for any density. 

3. w(n) -*7>0asn->oo: Condensation depends on the first order correction to w(n). For example , if 

w(n) =7(l + 6/n + (9(l/n 2 )) (13) 
then condensation occurs for large enough densities iff b > 2. 

Let us discuss case 3 further. First of all we note that the leading term of w implies that the process (a), see Fig.^ 
occurs with a rate proportional to nj_inj+i. This is reminiscent of the filament current @ which is proportional to 
the product of the filament density at two different locations. To discuss this case further, we will first consider 



(14) 



w(n) =7(1 + 



2 In general, the first order correction to 7 is not 0(l/n). A detailed study of these non-trivial cases can be found in l27l . 
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FIG. 2: Probability P(n) of finding in the steady state a site with n particles obtained from stochastic simulations. The different 
curves correspond to different models within the general class where the non-zero rates are specified by Eqs. JjJ-JSJ. Top: ZRP 
(u = v = 0), Middle: Exactly solvable case (all rates non-zero), and Bottom: v = 0. In all cases, the probability distribution 
shows first an algebraic decay and then a high probability for finding a site containing a high fraction of the particles. The 
non-vanishing parameters are in all cases r = 0.24, p = 0.6, q = 0.4 and w(n) = 1 + 2.5/n. The length of the system was 
L = 2000 and the number of particles were N = 6L. 



where the higher order terms are exactly 0. Then 

n _i 

fc=i 

which is asymptotically f(n) rj 7 _n n~ b . 

The distribution must respect the global conservation of particles, i.e., J2i n i — N. It is helpful to consider the 
system in the grand canonical ensemble, where the density p is controlled by the fugacity z. The grand partition 
function is 

F (*) = E (-)"/(") ( 16 ) 

which is well-defined for < z < z c , where z c is the radius of convergence of F(z). For the rate l|14f) one gets 
z c = 7 |27j - Knowing that the density of particles in the grand canonical ensemble is p{z) = (n) = z d] ? z , the 
condensation can be understood as follows: if p diverges for z — > z c then one can obtain any arbitrary density by 
tuning z. However, if p(z c ) < oo then the maximum achievable macroscopic density is just p c = p{z c ). Thus, in a 
system having density p > p c , the extra particles, in total (p — p c )L, must form the macroscopic condensate. For 
the rate H14[) it can be shown analytically |2jj that p c = is finite for b > 2 . For any other form of w(n) which 
asymptotically reduces to Eq. (|14fl . for example w(n) = 7(1 + b/(n + c)), b c is still 2 whereas p c is in general different 
from -r^K. When b > 2 and p> p Cl the probability of finding n particles at a site decays as n~ f ' with average density 
p c . The condensate occurs on top of this "critical fluid". See Fig. [21 for an example. 

Thus, b is the important parameter of the model which decides if the system can possibly phase seperate. Physically, 
in the ZRP, b can be interpreted as the escape rate of particles from the condensate j^]. Note that the mean rate at 
which particles hop from the "critical fluid" is given by (w) = 7 and thus particles escape from the condensate with 
a rate w(n) — (w) — jb/n. 

Comparing these results to the minimal model of active filament bundles we thus find a remarkable difference. 
While condensation depends in both cases on a critical density, in the hopping model the value of this critical density 
does not depend on the leading order term, but on the first order correction b. In contrast, the critical value of the 
minimal model depends on a which corresponds to the leading order term as will be shown below. 
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FIG. 3: Mean-field analysis of the stochastic model with rates specified in Eqs. 0-© with w(n) = 1 + b/n. The homogeneous 
state is stable for values of b below the blue line if v = and below the red line if v ^ 0. The parameter values are p = q = r = 1, 
7=1, and the limit L — *• oo is considered. 



In order to further investigate the last point we will now consider the case v = when there is no repulsion of 
particles. In this case an analytic solution is out of reach. Numerical solution reveals condensation for b > 2. As 
an example in Fig. we show the steady state distribution of particles for b = 2.5. It shows a clear peak at high 
particle numbers similar to the corresponding distributions for exactly solvable cases. For smaller particle numbers, 
the distributions show a power law with an exponent £ = 2.4 ± 0.05. The exponent is thus close to 6, which is the 
value obtained in the exactly solvable case. 

Some insight into the condensation behavior of the model can be gained by using a mean-field approximation. To 
this end we first write the dynamic equations for the expectation values Vi — (n^) for i = 1, . . . , L. These equations 
depend on expectation values of products of occupation numbers of different sites. The mean-field approximation then 
consists of approximating these by products of the expectation values such that the dynamic equations are closed. 
Following this procedure we find 



First, consider again the rates ©-© with w given by Eq. 1)14(1 . but now with v = and p = q = r = 1. It is easy 
to check that Ui = po = const, is a solution of ((T7jl . Because of the global conservation of particles Lp = N, thus 
po must be the density of the system. Stability of this homogeneous state v i = po can be checked by adding a small 
perturbation 8vi. Representing the density by a Fourier-series we find for the dynamic equations in linear order in 
the perturbation 



IV. ABSENCE OF PARTICLE REPULSION v = : MEAN-FIELD ANALYSIS 



dt 



dii(vi-i)vi-\ - d R {v i )v l - d L {vi)v t + d L (u l+1 )u l+1 

-u(l/i-2, Vi)Vi-2Vi + 2u(l/j_i, V i+1 )Vi-ll/ i+1 ~ u(vi, V l+2 )V l V i+ 2 



(17) 



d8v k 
dt 



(18) 



where 



e fe = 2 7 (cos(27rfc/L) - 1)(1 - 2 7 (6 + p Q ) cos(27rfc/L)) 



(19) 



and Srik are the Fourier coefficients of the perturbation, k — 1,...,L. Thus, the pertubation Sv^if) = e Cfct (5;// c (0) 
decays if > and grows in the opposite case. Since cos(27rfc/L) — 1) < for all k, modes can only become unstable 
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FIG. 4: Steady state density profiles of the mean-field equation l|17fl. Rates are given by 0-© with p = q — r — 1 and 
w(n) = 1 + b/n. Density profiles for a super- and a sub-critical value of b are shown. Parameters 7=1, po = 0.25, such that 
b c = 0.25, and b = 0.24 (blue line) and b = 0.26 (red line). 

if 2(6 + po)7 > I- In the limit L — > oo the instability occurs at b c = I/27 — po- The minimal model discussed in 
section III Al is equivalent to the 6=0 case of Eq. (|14f> . where the critical value of 7 is inversely proportional to the 
density pq. 

In the case when v is chosen according to Eq. © the critical value of 6 is given by 



where A = p — 1. Figure [3] describes the phase diagram in (6, po)-P ia ne for 7=1. Inclusion of the u-process shifts 
the critical line to larger values of 6, that is, it increases the stability region of homogeneous state. 

A few points resulting from the mean-field analysis merit attention. First, it is interesting that the condensation 
transition survives in the mean field approximation, which is not the case for ZRP, i.e. when v = = u. One can easily 
check that in this case = 2j(cos(2ir k/L) — 1) < and thus 6 C = 00. Second, although the mean-field equations of 
our model truly capture the transition, they do not reproduce the exact critical point 6 C = 2. Also, the density profile 
obtained in the mean-field approximation is different from the exact results, see Fig. 0J In the stochastic system, the 
condensate appears on the top of a critical fluid, whereas the mean-field equations hide the fluctuations coming from 
the fluid and all the matter is found in the condensate. It is also interesting to note that the mean-field analysis would 
not change if d is taken to be a constant d = 7 instead d = 7(1 + 6/71), because first line of Eq. I|17|) gives identical 
results for both cases. This interesting case d = 7 is further explored in next section. 

V. CONNECTION TO PHENOMENOLOGICAL DESCRIPTIONS OF ACTIVE FILAMENT BUNDLES 

Having discussed the driven diffusive system introduced in Section [H] for a specific choice of the rates, we now want 
to return to the general case and to clarify the connection of the stochastic models to the minimal model of active 
filament bundles which motivated our analysis. This will be done through the intermediate of phenomenological 
equations governing the dynamics of such bundles on large length scales. 

Phenomenological theories of active gels - of which active bundles are a special example - are based on the observation 
that filament currents are generated by gradients in the stress 0, El El The system's stress can be expressed in 
terms of the state variables, that is the dynamic fields describing the system. If all filaments are aligned in a bundle 
and point all into the same direction, the stress tensor reduces to a scalar a that depends on the filament density only. 
Furthermore, if any elastic response of the system is neglected, then the relation between the stress and the filament 
current j can be written as 



(A 2 + 1/2) + + 1/2) 2 - 2A3 + (1 - 7 - 1 M 2 , 



(20) 



j = T]d x a 



(21) 
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where r\ is an effective friction coefficient [l^. ^ n the phenomenologicai descriptions, the stress a is then systematicaliy 
expanded in terms of the filament density and its gradients with respect to some homogeneous reference state. 

We can connect the hopping model to the phenomenologicai dynamic equation by considering a coarse grained 
version of the discrete mean-field equations of the hopping model. The coarse-graining procedure amounts to assuming 
that the occupation numbers of the lattice sites vary only weakly with the site index i. On large length scales the 
discrete index can then be replaced by a continuous variable x and the occupation numbers Vi can be replaced by a 
density c such that c(x = id) = Vi/a. Here, a is the lattice spacing. Replacing v± in the discrete mean-field equations 
and then expanding terms like c(x ± a) into a Taylor-series c(x ± a) — c(x) ± c'(x)a + c" ' (x)a 2 ± . . ., which is truncated 
at some order of a, we arrive at a partial differential equation for the time-evolution of the density. 

We now apply this procedure to the mean-field equations of the stochastic model (i.e, Eq. I17f) . where we focus for 
simplicity on cases when only particles of one kind are present and when = du = d. Restricting the expansion to 
terms of at most 4th order in a we find 



d f c = 



-dl 



h(c) + -d 2 x h(c) 
7 

9(c> c) + —dlg(c, c) - 2did 2 g(c, c)(d x c) 2 



(22) 



where without loss of generality we have set a = 1 in the final expression. In this equation, the hopping rates 
are contained in the functions h and g. Explicitly, h(c) = d(c)c + v(c)c and g{c,c) = u(c,c)c 2 . Furthermore, 
d t g(c,c) = d Cl g{ci > C2) |d =c 2 =c The dynamic equation for the filament density derived from the discrete mean-field 
equations of the hopping model are thus of the same form as the corresponding equation for the filament density in 
the phenomenologicai descriptions. From this we can deduce the important result that the sum of the hopping rates 
with appropriate sign multiplied by the filament density equals the tension <jq m the homogeneous state 



(To = u(cq, c )cq - (d(c ) + v(cq))c . 



(23) 



This relation provides a simple link between the kinetic hopping rates, which are microscopic quantities, and the 
stress generated in the bundle, which is a macroscopic quantity. It should be noted, however, that this relation is in 
general one way. While for every microscopic model the stress can be determined, the opposite is usually not true. 

As an example consider the rates introduced in Section llVI namely d(c) = D, v(c) = 0, and g(c, c) = 7 2 (c + b) 2 . In 
this case the dynamic equation reads 



d t c= {D -2b 1 2 )d 2 x c- 



D - Ubj' 
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di 



l 2 dl 



G 



(24) 



This equation can be compared to the expression obtained in the phenomenologicai approach developed in Ref. |15| . 
This comparison shows that all terms appearing in the general expansion are generated by the stochastic model. It 
can therefore be used to study the effects of fluctuations on the dynamics of a generic active filament bundle. Note, 
that in Ref. [15] the term cx d x c was neglected for reasons of simplicity. 

A linear stability analysis of Eq. 124f) yields a critical value b c of the parameter b for which the homogeneous 
distribution loses its stability. This corresponds to condensation in the stochastic model. For the critical value we find 
b c = D/2j 2 — Co. This is the same value as obtained in Section Hvl for L — » 00 in the mean- field limit of the discrete 
hopping model, where D = 7. The corresponding tension created in the homogenous state at the critical point can 
be expressed as 00 = 7 2 (co + b)(b + c — 2). An instability of the homogenous state leads to a state where filaments 
condense, i.e., accumulate at one point and can therefore be interpreted as the point where the bundle ruptures. The 
value of the tension at the critical point is thus the maximal tension that can be generated by a filament bundle before 
it tears apart. 



VI. DISCUSSIONS AND CONCLUSIONS 



In summary, we have introduced a new class of hopping models on one-dimensional lattices. Motivated by the 
dynamics of filaments connected by molecular motors, this class is defined by processes in which two particles hop 
simultaneously; either from adjacent sites to a site in between or away from one site to both adjacent sites. This is, 
in particular, different from the ZRP where such processes occur in an interval At with a probability tx (At) 2 and 
are therefore neglected in comparison to processes that occur with a probability cx At. As a consequence of the two- 
particle processes the mean-field approximation shows like the stochastic system a condensation transition. This does 
not hold for the mean-field analysis of the ZRP. Despite these important differences, for rates that have a factorized 
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steady state probability distribution, the steady state distributions of the ZRP and of the class of systems introduced 
here are identical. This similarity allowed us to transfer results from the ZRP on the conditions of condensation to 
the new model class. In particular, we found that there is a critical density for condensation, if the rates approach 
for large occupation numbers a constant different from zero. 

Coarse-graining of the mean-field equations allowed us to connect the kinetic hopping rates in the stochastic model 
to the stress generated in the bundle. The microscopic rates are thus linked to macroscopic physical quantities. Apart 
from giving the rates a physical meaning, this connection could be used to obtain the strength of the noise term in 
Langevin equations for the time-evolution of the filament density. As active filament bundles are inherently out of 
thermodynamic equilibrium there is in general no fluctuation-dissipation theorem that would allow to determine this 
strength. The coarse-graining also points to limitations of continuum descriptions of active filament bundles. Indeed 
we find that the criterion for condensation in the exactly solvable case involves only the first order correction of the 
rates, while the critical value in the mean-field equations depends also on the zeroth order parameter. 

Taken these points together, the discrete hopping models introduced here are seen to be a useful tool for analyzing 
the effects of fluctuations on dynamics of active filament bundles - in spite of their simplicity. In the future it will 
be in particular interesting to study the case of two different particles present on the lattice. In the corresponding 
microscopic and phenomenological descriptions, traveling waves were found in this case. 
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